Lorentz shear modulus of fractional quantum Hall states 
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We show that the Lorentz shear modulus of macroscopically homogeneous electronic states in 
the lowest Landau level is proportional to the bulk modulus of an equivalent system of interacting 
classical particles in the thermodynamic limit. Making use of this correspondence we calculate the 
Lorentz shear modulus of Laughlin's fractional quantum Hall states at filling factor v — 1/m (m an 
odd integer) and find that it is equal to ±h~mn/A, where n is the density of particles and the sign 
depends on the direction of magnetic field. This is in agreement with the recent result obtained by 
Read in arXiv:0805.2507 and corrects our previous result published in Phys. Rev. B 76, 161305 
(R) (2007). 
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I. INTRODUCTION 

The response of many-body systems to external fields 
can be universally formulated in terms of hydrodynami- 
cal equations of motion, which follow from the local con- 
servation laws of the number of particles and momentum. 
These equations completely describe the dynamics of the 
particle density n and the particle current density j under 
the combined action of external and internal (or stress) 
forces. The most important ingredient of this formula- 
tion is the stress tensor P^, which determines the stress 
force density entering the momentum conservation law: 
^stress _ —Q k p ik j n general the stress tensor for ev- 
ery particular state of matter is a universal functional 
of the current density^ (or, equivalently, of the deforma- 
tion tensor^). In fact, the form of Pik[j) encodes all the 
dynamical information about a given many-body state. 

Unfortunately it is an extremely difficult, if not hope- 
less task to determine the general form of the stress ten- 
sor as a functional of the current. As usual the prob- 
lem simplifies in the linear response regime as the stress 
tensor becomes a linear functional of the strain tensor 
u ik = ^(diUk + dkUi), where u(r, ui) is the displacement 
vector field. The coefficients of this linear functional form 
the rank-4 tensor of elasticity Qijki{w), which is, in gen- 
eral, a function of frequency. This tensor describes a 
stress response to a time dependent deformation of the 
system. In a macroscopically isotropic two-dimensional 
(2D) system subjected to a perpendicular magnetic field 
the rank-4 tensor Qijki contains only three independent 
components and can be parametrized by three dynamic 
"elastic moduli" K(w), n{w), and A(w)£MiI. The mod- 
uli K (oj) and n{uS) describe the response to a local change 
in volume and to a local volume-preserving (shear) defor- 
mation, respectively. Therefore they correspond to the 
standard bulk and shear moduli of classical elasticity the- 



ory4 The third modulus, A(w), appears only in the pres- 
ence of the magnetic field and plays an important role 
in the dynamics of 2D electrons at high magnetic field. 
This modulus controls the magnitude of a stress propor- 
tional to the rate of volume-preserving deformations, and 
in this respect is similar to a viscous stress. However, the 
corresponding force always acts in a direction perpen- 
dicular to the stream lines and is purely nondissipative, 
like the usual Lorentz force. To underline its nondissi- 
pative character, in Refs. HH0 we named the modulus 
A(w) "Lorentz shear modulus". The same quantity is 
also known in the literature as "asymmetric viscosity"—, 
or "Hall viscosity"—. 

The appearance of a "nondissipative viscosity" in the 
presence of a magnetic field is well known in plasma 
physics. The calculation of the corresponding kinetic co- 
efficient Ao = — lim w _>o A(cj) for a classical plasma can 
be found, for example, in Ref. I 111 . The result is 



A = ± 



hn k B T 

2 hhJ r 



(1) 



where T is the temperature and u> c is the cyclotron fre- 
quency (the sign is determined by the direction of the 
magnetic field). 

In this paper we focus on the microscopic calculation 
of Ao for an extended, macroscopically homogeneous 2D 
electronic system in the quantum Hall regime, i.e. when 
all the electrons reside in the lowest Landau level. For 
a 2D noninteracting electron liquid in a completely filled 
Landau level, this calculation was first done by Avron et 
al. in Ref. 9. These authors observed that the value of 
Ao is proportional to a Berry curvature related to adia- 
batic changes of geometry, and used this fact to explicitly 
calculate the modulus. The result was A = ±frn/4ji2 
which, incidentally, agrees with the classical formula ([I]) 
if one replaces the thermal kinetic energy per particle 



2 



ksT with the kinetic energy per particle of the lowest 
Landau level Hu) c /2. 

In Ref. one of us derived hydrodynamics equations 
for quantum Hall states at fractional filling factor us- 
ing a fermionic Chern-Simons theor y 13 ' 14 i 15 i 16 at the 
RPA level. Within this approximation Ao exactly co- 
incides with the noninteracting result of Ref. |9|, as re- 
ported above, which is not surprising since the compos- 
ite fermions do not interact in RPA. In a recent paper— 
we applied a formally exact linear response theory to the 
calculation of the Lorentz shear modulus for strongly cor- 
related quantum Hall states at fractional filling factors 
(Laughlin states). We showed that A(w) can be expressed 
in terms of a particular stress-stress correlation function 
which, in the limit lu — > 0, reduces to the Berry curva- 
ture expression of Avron et ai. After publication of our 
work 7 the problem of calculating Ao was reconsidered by 
ReacU* who noticed an error at the very end of our cal- 
culation for Laughlin states. Making use of the Laugh- 
lin plasma analogy 17 Read correctly calculated the Hall 
viscosity (Lorentz shear modulus) for Laughlin states at 
v = 1/m (m is an odd integer) and also for states with 
trial wave functions in the form of conformal blocks of a 
conformal field theory. For Laughlin states, Read's result 
is 



whereas in our paper we had incorrectly obtained Ao = 

I Hn 
4 ' 

Read's paper is difficult and contains much more than 
just the calculation of the Lorentz shear modulus. We 
have found that it is possible to give a more elementary 
derivation of Read's result for filling factors v = 1/m. 
The crucial step in the derivation is the recognition that 
the calculation of the Lorentz shear modulus for a macro- 
scopically homogeneous state of electrons in the lowest 
Landau level can be mapped to the calculation of the bulk 
modulus of an equivalent system of interacting classical 
particles. While in general the many-body interactions 
in this equivalent classical system are prohibitively com- 
plicated, they simplify dramatically for Laughlin's quan- 
tum Hall states, where one recovers the well-known clas- 
sical plasma analogy. This allows a simple calculation of 
the Lorentz shear modulus. The quantum-classical cor- 
respondence can also be used "in reverse" . Namely, from 
the knowledge of the shear modulus associated with a 
certain wave function in the lowest Landau level one can 
in principle obtain the bulk modulus of the equivalent 
classical system, even if the latter is very complicated. 

The structure of this paper is as follows. In Sec. IIA 
we introduce the formal definition of the Lorentz shear 
modulus an discuss its physical significance. Here we also 
address the delicate question of the correct boundary con- 
ditions for the trial many-body wave function to be used 
in the calculation of the Berry curvature in the thermo- 
dynamic limit. The calculation of the Lorentz shear mod- 
ulus for the Laughlin states is presented in Sec IIB. In 



Sec. IIC we generalize our procedure to other macroscop- 
ically homogeneous states in the lowest Landau level. In 
Sec. Ill we summarize our main results and discuss how 
the new value of Ao affects our previous analysis of the 
collective modes of the electron liquid in the fractional 
quantum Hall regime. 



II. LORENTZ SHEAR MODULUS IN A 2D 
MAGNETIZED ELECTRON GAS 

A. Definitions and physical significance of the 
Lorentz shear modulus 

Let us consider a 2D electron gas confined to the (x, y) 
plane and subjected to a perpendicular magnetic field 
B = Bz. In the linear response regime and in the long 
wavelength limit the exact stress tensor takes the 
following "elastic" formal 7 - 

Pij(r,ui) = -Qijki{u)u k i{r,u)), (3) 

where um = \{dkUi + diUk) is the strain tensor, and u 
is the displacement vector defined in the standard way, 
namely <9 t u = j/n = v is the velocity field of the electron 
liquid. The frequency dependent coefficients Qijki{u) m 
the linear functional of Eq. |3|) form the dynamic tensor 
of elasticity. The general structure of this tensor is essen- 
tially fixed by symmetry. In particular, for a macroscopi- 
cally isotropic state the rank-4 tensor Qijki is completely 
determined by only three independent "elastic moduli" 

Qijki(w) = K{Lo)5i 3 8ki + fi(u))(6ik6ji + SuSjk - SijSki) 
+ iw y~ (sikSji + EjhSu + suSjk + £ji6ik)(4) 

Inserting Eq. © into Eq. ^ we get the folowing repre- 
sentation for the stress tensor for any isotropic state of a 
2D magnetized electron gas 

P%j = - KdijUkk - M( 2 «ij - SijUkk) 

+ A(e ik v k j +SjkVki), (5) 

where Vij = —iuiUij = \{diVj + djVi) is the rate of change 
of the strain. The first two terms in Eq. ([5]) describe the 
stress arising from a local change in volume (ukk = Vu) 
and from a volume-preserving (traceless, or shear) de- 
formation (2uij — S^Ukk), respectively. Hence the cor- 
responding coefficients K(u>) and fJ>(u>) have the same 
meaning as the standard bulk and shear moduli of a ho- 
mogeneous elastic medium^ The third term in Eq. ([5]) is 
proportional to the rate of change of the strain, i.e., to 
velocity gradients, which looks somewhat similar to a vis- 
cous stress. However, for real A this term is time-reversal 
invariant and, therefore, does not cause dissipation. To 
get a more intuitive picture we consider a "shear flow" 
with the velocity distribution of the form v = (v x (y),0) 
(for example, two oppositely directed streamlines). In 
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p xx= A ^,"Ay + dy) 



P yy =Ao3,v I (.v) 



i 



T 



v x (y+dy) 



Pyy= Ao3,V r (v) 



v x (y-dy) 

P xx = A r) \ i : . ,h > 




FIG. 1: Lorentz shear forces acting on a small rectangular 
element of the electron fluid according to Eq. ©. 



this simple case the A-part of the stress tensor has only 
two nonzero elements, 



P 1 : 



-PS; 



Ad y v x (y). 



(6) 



This distribution of stresses corresponds to forces, ex- 
erted on the faces of an infinitesimal rectangle, which 
squeeze or repel (depending on the sign of A) two op- 
posite stream lines. The important point is that these 
forces are always perpendicular to the local direction of 
the streamlines, as shown in Fig. 1. Hence the effect 
of the third term in Eq. ([5]) is analogous to that of the 
Lorentz (Ampere) force acting between two currents, and 
causes no dissipation. Having in mind this physical pic- 
ture we prefer to call the coefficient A(w) the "Lorentz 
shear modulus " in contrast to more formal terms such as 



"asymmetric viscosity 
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"Hall viscosity 
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Linear response approach to the microscopic 
calculation of the Lorentz shear modulus 



Let us now address the problem of the microscopic 
calculation of the Lorentz shear modulus in the limit of 
zero frequency. In our recent paper— we expressed Ao = 
— lim^^o A(w) in terms of the stress-stress correlation 
function in the following manner: 



Ao = — lim lim —lm((P x 



(7) 



The stress tensor operator P^ is formally defined as fol- 
lows (see, e. g., Ref. (3) 



Pij(r, t) 



5H[ 9l 
Sgij(r, t) 



J gij=Sij 



where H[gij\ is the Hamiltonian of the system in a "de- 
formed" space with metric <?jj(r,i). The subscript q in 
Eq. ([7]) means that we are actually considering the corre- 
lation function for the Fourier component of Pij{v,t) at 
wave vector q, where q tends to zero before lo. It is con- 
venient to parametrize the metric tensor as follow o 7 ' 9 ' 18 



J 

T2 



1 Ti 



n 



(8) 



FIG. 2: Deformation of the Euclidean plane corresponding 
to the metric of Eq. ((8|. The square is transformed into a 
parallelogram. The slope of the oblique side coincides with 
the direction of r in the complex plane. We have chosen 

J = 1, Tl = T 2 < 1. 



where r = t\ + it-i- A picture of the deformation of the 
Euclidean plane corresponding to this choice of metrics 
is shown in Fig. 2. Making use of the Lehmann repre- 
sentation-^ for the stress-stress correlation function and 
performing the standard manipulations one can trans- 
form Eq. (JT)) to the following Berry curvature form 



A °=^ im w 



3*0 
<9t 2 



(9) 



where L 2 is the area of the system, ^>q is the ground state 
wave function in a homogeneously deformed space with 
a constant metric of Eq. ((8|), and the r-derivatives are 
calculated at t\ = 0, T2 — 1, and J = 1. 

Equation © was first derived in Ref|Ej from the adi- 
abatic response theory. Our linear response derivation 
shows that there is a delicate point in the identification 
of the Berry curvature [the right hand side in Eq. 
with the physical Lorentz shear modulus. Physically such 
a modulus describes a stress response to a deformation 
whose wavelength is much larger than any internal scale 
of the system (e.g. the interparticle distance, the correla- 
tion length, etc.), but still much smaller than the size of 
the sample L. At the level of the linear response formula, 
Eq. ([7]) , this means that the thermodynamic limit should 
be performed before the limit q — > 0. This guarantees 
that the calculated elastic modulus does not depend on 
the geometry of the sample (provided the sample is suf- 
ficiently large), i.e. it is a bulk property of the "ma- 
terial", as it should. However, in the derivation of the 
Berry curvature formula the order of limits was explic- 
itly interchanged. The right hand side of Eq. © con- 
tains the ground state wave function of a finite system 
with a finite number of particles N. Physically this de- 
scribes the response to a homogeneous deformation of the 
whole sample which in general may depend on the sample 
geometry^ Nonetheless it is still possible to obtain the 
correct physical result even with the interchanged order 
of limits. The key is to use boundary conditions for the 
wave function which are compatible with the symme- 
try of the physical plane wave perturbation (the defor- 
mation in our particular case). For, in this way we ensure 
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that the q — > limit of the deformed state goes smoothly 
to the homogeneous deformation of the ground-state in 
the same geometry. In practice this simply means that 
the correct result is guaranteed if one studies a sample of 
rectangular shape, which satisfies periodic boundary con- 
ditions either in two directions (torus) or in one direction 
(cylinder). On the other hand a sample of circular shape 
would not be compatible with the symmetry of the plane 
wave, because the presence of a finite wave vector, no 
matter how small, always breaks rotational symmetry 22 

The choice of the torus geometry in the present context 
may also cause some technical problems which require a 
special care. On a torus the ground state for gapped 
quantum Hall states is degenerate. This nicely demon- 
strates the topological nature of these states, but at the 
same time introduces an artificial degeneracy which is not 
present in any real physical sample. Clearly, this degen- 
eracy should not affect the physical elastic moduli. To 
guarantee that this is indeed the case we put our system 
of N electrons on a cylinder with a circumference L in 
such a way that it occupies the area L 2 = Ns(2ir£ 2 ), 
where N$ — N/v is the number of flux quanta (one 
flux quantum = hc/e), and £ — ^/Hc/eB is the mag- 
netic length. The thermodynamic limit which we take at 
the end of calculations corresponds to N —> oo at fixed 
filling factor v = N/N s . 

For the Hamiltonian defined on a space with the metric 
of Eq. ((8]) the TV-body wave function which lies entirely 
in the lowest Landau level and satisfies periodic bound- 
ary condition in the x-direction can be written in the 
following general form 



normalization integral to the following form 

L 



Z(r 2 ) 



- 2 - n n Ux k 

k=i i "L 



dyk 



o 

J2-7TIZX/L 



2mz N /L 



— J,, 

r, 2^3 = 1 t,( ; 



(12) 



where Zk — Xk + iyk- Hence the final formula for the 
Lorentz shear modulus in the state \&n takes the form 



An = 



2L 2 



d 2 

M lnZ{T2) 



(13) 



7-2 = 1 



with Z(t 2 ) defined after Eq. (fl2f . In the next subsec- 
tions we calculate this integral using the Laughlin classi- 
cal plasma analogy. 



C. Calculation of the Lorentz shear modulus for 
Laughlin states 

The cylindrical generalization of the Laughlin trial 
function at v — 1/m contains the analytic factor / of 
the following formal 



,**)=n( 



^2-KiZj/L ^2iriz k /L 



(14) 



]<k 



Inserting this equation into Eq. (|12p we represent the 
normalization integral in a form of a partion function of 
TV classical particles, 



y (rx,...,r N )=Z 2 (T 1 ,T 2 )f(r) 1 



N 

IT 



T T Vi 



(10) 



where Z 2 is the normalization factor, 



exp(2Tri^-ip^), and /(771, . . . , ?yjv) is an analytic func- 
tion of its arguments. A crucial observation is that \tn, 
Eq. (|10j) . apart from the normalization factor, is an an- 
alytic function of the complex variable r. This enables 
us to express the Berry curvature solely in terms of the 
normalization factoriii& 



21m 



a* 



dr 2 



dr 2 



InZ. 



(11) 



N 

Z(t 2 ) = T2 N f[ J dx k J 



dy k e 



-;ff(n,.,r») 



(15) 



at the "temperature" v = 1/m and with the following 
energy 

2irN N 1 
W(t u ■■■,r N ) = -j 2 J2y^ + 2^ V(r 3 ,r k ), (16) 

« A 1 _L J.. 



i=i 



where we used the identity 1/m = N/N s = 2irN£ 2 /L 2 , 
and introduced the notation 



V (r,r'l 



hi 



(17) 



Hence the problem of microscopic calculation of the 
Lorentz shear modulus reduces to the calculation of the 
normalization factor. For our particular choice of a quan- 
tum Hall system on a cylinder the problem simplifies even 
further since the normalization factor depends only on r 2 . 
Indeed, by shifting all the a;- variables, (xj + Tij/j) — > Xj, 



and rescaling all y- variables, r 2 yj 



y.i 



we reduce the 



The first term in Eq. (|16p has a clear physical interpre- 
tation: it is the Coulomb potential of a homogeneously 
distributed (on the cylindrical surface) positive charge 
with density p = N / (t 2 L 2 ). Note that after the rescaling 
of the j/-coordinate the size of the system along the axis 
of the cylinder becomes L y = t 2 L: the classical system 
with the probability distribution ~ exp[— mH/({rj})] oc- 
cupy the region with < y < t 2 L. Hence the quantity 
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N / (T2L 2 ) entering the first term in Eq. (fT6f exactly coin- 
sides with the physical density of this classical system, in 
perfect agreement with the idea of the Laughlin classical 
plasma analogy. 17 However, at variance with the original 
circularly symmetric Laughlin wave function, the second 
term in Eq. IT6|) does not look like the interaction energy 
of a Coulomb plasma. To show that this is nonetheless 
the case, we express the two-point "potential" V^r, r'), 
Eq. (fT7|) . in terms of the physical Coulomb interaction 
V(r — r') between two point particles on a cylinder, (see 
Appendix A) 



V(r-r') = -2tt 



\y-y 



-In 



l-e 2 



-tc' + i 



(18) 




FIG. 3: Distribution of the electron density in the y direction. 
Notice the two small edge regions, whose width is proportional 
to the mean interparticle distance i ~ -7= . 



The identity of Eq. (1A3|) allows us to relate V(r, r') to 
V^r-r') 

V(t,t') = 2n^-^- + V(r-r'). (19) 



Substituting this equation in to Eq. (JT6J) we transform 
the energy ^({rj}) to the folowing form 



W({r 3 }) 



-r 2 ^(/v- 



l) 2 



2ttN 



N 



1 



t 2 L 2 ^ 



-rfe) 



(20) 

where Xf. = and y^ = yf. + \t 2 L{\ — jj). Hence the 
main effect of the difference between V(r, r') and V(r— r') 
is a shift of the y-coordinate by a half of the system size 
~ T2L/2. In what follows we remove this shift by moving 
the origin to the center of the slab: y% — > y^. 

Thus the energy W({rj}) exactly corresponds to the 
energy of a system of N classical charges on an infinite 
homogeneously charged cylinder. To proceed further in 
the calculation of the partition function of this system 
we separate a "neutralizing" part of the background with 
density 



n b (r) 



N 



L 
T2 2 



\y\ 



(21) 



and total charge J drn^r) = N. The next step is to sep- 
arate the energy of interaction of the particles with rij,(r) 
from the total energy of Eq. (|2H|) . Using Eqs. (|A4[) and 
(|A5|) of Appendix A we find the following representation 
for Wdrj}) 

f2N 2 1\ 
W = -T 2 ir i — -N+-J+ E NJ {{rj}) + E s ({rj}). 

The second term, E^j({rj}), in this equation is the en- 
ergy of a neutral slab of jellium, 

E ^ = \J dvdv'V{v - r')[A ( 5(r)A ( 5(r') - <5(r - r')p(r)] 

(23) 



where p(r) is the microscopic particle density, and Ap(r) 
is the microscopic charge density of jellium: 



N 



p(r)=J2s(r-r J ), Ap(r) = p(r) - n 6 (r). (24) 



The third term in Eq. 
(edge) contribution 



22|) corresponds to the surface 



^ = -^2.(1^1- — 



bl- — 



(25) 



which is the potential energy of particles outside the 
"neutralizing" part of the background (see Fig. 3). 

Eq. (|22|) is extremely useful for the analysis of the rel- 
evant limit N — > 00 of the partition function. It it phys- 
ically obvious that the contribution, Eq. (|25|) . should be 
irrelevant in the thermodynamic limit as it involves a 
small fraction of particles in the edge region. Indeed, 
from Fig. 3 we see that the contribution to E$ comes 
from a region whose width is of the order of the mean 
interparticle distance I ~ -4= = 77^ (since it is propor- 
tional to the magnetic length we use the same notation) . 
Then a simple, order of magnitude estimate of E$ is 



E S 



N 



(L£n)£ 2 



N . 



(26) 



where Lin ~ v N is the average number of particles in 
the edge region. This shows that in the limit N — > 00 
Es grows as V^/V and is negligible in comparison with 
the first and the second terms in Eq. (|22|) . which are 
proportional to N 2 and N respectively. Inserting the 
energy of Eq. (|2"2"|) into Eq. (fT5|) and neglecting Es, we 
rewrite the normalization integral as follows 



Z(T2. 



_ T 2N rnT 2 n 
— lj 6 



2 .v 
3 



2b ' 



(27) 



where .Znj = (exp[— m.E?Nj({ij})]) is the partition func- 
tion of the slab of jellium, and Z = (t 2 L 2 ) n is the par- 
tition function of TV noninteracting particles in the slab. 
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Hence Eq. 
form 



27]) can be equivalently represented in the 



Z(t 2 ) = L ™ e mT 'H^- N+ ^ e - mF "^\ (28) 

where Fmt = — \u(Zn.i /Zo) is the interaction free en- 
ergy of the classical jellium model. The final step is to 
substitute Eq. (|28|) into the expression of Eq. (fl~3|) for 
the Lorentz shear modulus. The first, non-extensive ex- 
ponent (~ A^ 2 ) in Eq. (f28|) does not contribute to Ao 
since it is linear in r 2 . Therefore Eq. (fTB"]) reduces to the 
form 



Ao = - 



hm 



d 2 F int (r 2 ) 

— T 2 







T-2 



L 2 



(29) 



Now we can easily perform the thermodynamic limit. 
In this limit, the free energy density of a neutral jel- 
lium depends only on the density of particles, which is 
equal to the background charge density p of the jellium: 
F N j/(t 2 L 2 ) = / N j(p), where p = n/r 2 and n — N/L 2 is 
the density of electrons in the original quantum Hall sys- 
tem. Thus in the thermodynamic limit Eq. (|29p simplifies 
as follows 



A 



hm 
~2 



Ot: 



Tlfii 



(30) 



where /i n t(p) is the interaction free energy density of a 
2D classical jellium model with logarithmic interaction 
between the particles. The key observation is that t 2 
enters Eq. (|30[) as an effective volume (the area in 2D). 
Hence the second derivative of the energy with respect to 
r 2 in Eq. (|30|) is nothing but the isothermal bulk modulus: 



A 



hm 



d_ 

dr 2 



Pn 



hm 



-Ku 



(31) 



T 2 =l 



where P m t{p), and Ki nt (p) are the interaction contribu- 
tions to the pressure and bulk modulus, respectively. 
Equation (|31[) is the main result of the present paper. 
It shows that the calculation of the Lorentz shear modu- 
lus in a quantum Hall system reduces to the calculation 
of the bulk modulus of an equivalent classical system. In 
the next subsection we will demonstrate that this result 
holds more generally for any macroscopically homoge- 
neous quantum Hall state. 

For the Laughlin states considered in this subsection 
we need the bulk modulus K ln t(p) of a 2D Coulomb 
plasma with logarithmic interaction, which can be cal- 
culated exactly^ The easiest way is to use the virial 
representation for the interaction pressure 

P ^ = -JS J drdr'(r-r') d ^~^\ (r)p(r%(r,r')-l] 

(32) 

where S is the area of the system and g(r, r') is the 
pair correlation function. Inserting V(r) = — In |r| 2 into 
Eq. (|32p and using the sum rule 

J dr'p(r')[g(r,r') -!]=-!, 



we obtain the interaction pressure and the corresponding 
bulk modulus in the form 

Pint{p) = Kint{p) =-\p- 

This implies the following final result for the Lorenz shear 
modulus of the Laughlin quantum Hall liquids 



An 



hnm 



irh 
2P' 



(33) 



This equation reproduces the recent result of Read 10 . In 
the case of a full Landau level (m — 1) it recovers the for- 
mula for the "asymmetric viscosity" obtained by Avron, 
et. al. in Ref^. It is interesting to note a surprising fea- 
ture of Eq. (j33|) : because the density n is proportional to 
1 jm the Lorentz shear modulus for Laughlin states does 
not depend on the filling factor. 



D. Generalization to arbitrary macroscopically 
homogeneous quantum Hall states 

The result of Eq. (|3"Tj) can be straightforwardly gen- 
eralized to other macroscopically homogeneous states in 
the lowest Landau level. These are states in which the 
macroscopic particle density is uniform. For example 
a Wigner crystal, while microscopically inhomogeneous, 
can be considered homogeneous on a macroscopic scale: 
the average density is uniform. In fact, most states of 
interest in condensed matter physics have this property. 

In general the normalization integral for any wave func- 
tion ^o({ r j}) that satisfies cylindrical boundary condi- 
tions in the lowest Landau level is representable in the 
form of Eq. (fl~5| , where 



W({rj}) 



2-kN 



N 



T-2 



L 2 



$>? + &({r,.}), 



(34) 



3=1 



and U({fj}) = — v In \f({Tj})\ 2 is the interaction en- 
ergy of the equivalent classical system [here /({r^}) is 
the analytic factor in the many body wave function, 
Eq. ([TO])]. The interpretation of the first term in Eq. ([34]) 
is exactly the same as for the Laughlin states: it is 
the energy of classical particles in the electrostatic field 
of a homogeneous charged background with the density 
p = N/(t 2 L 2 ). Depending on the form of the wave func- 
tion ^Oj the second term can contain both two-particle 
and multi-particle interactions. In general the form of 
these interactions is extremely complicated: they are 
neither translationally nor rotationally invariant. How- 
ever, if the physical density distribution of the quan- 
tum state under consideration is macroscopically homo- 
geneous with a mean density n = N/L 2 , then the in- 
teraction energy U({rj}) of the equivalent classical sys- 
tem must contain a long-range Coulomb contribution of 
the form J2j<k ^( r ii r fc)- The presence of this pairwise 
logarithmic interaction is mandatory to compensate the 
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background potential and protect the macroscopic homo- 
geneity of the classical gas. Hence for macroscopically 
homogeneous states in the lowest Landau level the en- 
ergy of the equivalent classical plasma should take the 
form 



W({ ri }) 



2ttN 



N 



•E 

j<k 



Vir^r^ + U^}), (35) 



where the last term U sr ({rj}) can contain only short 
range (possibly multiparticle, anisotropic, etc.) interac- 
tions. Since short range interactions do not spoil the 
extensive character of the free energy, the chain of ar- 
guments that led us from Eq. (| 1 5[) to Eq. ([30]) remains 
valid for any macroscopically homogeneous quantum Hall 
state. Therefore the Lorentz shear modulus for such 
states can be calculated as follows 



A 



h 

2^ 



\ d p 




0^2 lnt 







T 2 =l 



2v 



K int (n), (36) 



where Pmt(p), and i^int(p) are the interaction pressure 
and the isothermal bulk modulus of the corresponding 
classical plasma. 

As a simple example consider the wave function pro- 
posed by Wexler and Ciftja3£ for the nematic liquid crys- 
tal state at v = 1/3. The cylindrical generalization of this 
wave function contains an analytic factor / the form 



f(zi,...,z N ) = II (' 

(• 



27ziZj/L 2izizk / 



2iTi{zj —a/2)/L 



j<k 

2Tri(z k +a/2)/L 



' e 2TTi(z j +a/2)/L _ e 2m(z k -a/2)/L 



(37) 



where a is a complex number, which we set equal to 
ia, with a real and positive. With this choice the sys- 
tem remains invariant under rotations about the axis of 
the cylinder. The parameter a is the microscopic dis- 
tance (of the order of the magnetic length) by which two 
of the three zeroes of the v = 1/3 Laughlin wave func- 
tion are displaced from their "regular" position on top of 
the particle. The Laughlin wavefunction is recovered by 
setting a = 0. It is now straightforward to verify that 
the classical system that is equivalent to this wave func- 
tion contains, in addition to the interactions of Eq. (|16[) . 
a two-body short-range interaction of the form (in the 
physically relevant 2D regime, \z% — Zj\ <C L) 



Mfe}) = -- 



2 E 



In 



1 



(38) 



This interaction is strongly anisotropic and exhibits log- 
arithmic singularities at interparticle distance (inde- 
pendent of direction) and a (along the y axis). It is 
clearly very difficult to calculate the bulk modulus of this 
system from classical statistical mechanics. However, if 
we could, by some independent method, calculate the 
Lorentz shear modulus of the wave function (|37p. then 
the problem would be solved. 



III. CONCLUSION 

In this paper we have shown that the calculation of 
the Lorentz shear modulus in a macroscopically homoge- 
neous electronic system in the lowest Landau level can 
be mapped to the calculation of the bulk modulus of an 
equivalent classical system. Application of this approach 
to Laughlin's fractional quantum Hall states gives a value 
of the modulus independent of filling factor in agreement 
with the result of Ref. llCJ and in contrast with what we 
had erroneously found in Ref. 0. In that paper, the in- 
correct value of the Lorentz shear modulus was used in 
the calculation of the collective modes of the fractional 
quantum Hall liquid, and the resulting dispersions were 
found to be in good agreement with the ones obtained in 
Ref. 1 by a more phenomenological approach. 

A serious problem arises when we repeat the calcu- 
lation of the collective mode dispersions according to 
Eq. (21) of Ref. 0, but using the correct value of Ao 
for the v = 1/3 Laughlin state. We get an instability 
in which the upper collective mode has negative oscilla- 
tor strength. This means that the simplest single-mode 
ansatz for A(u>), Eq. (20) of Ref. fails. Most likely this 
indicates that in reality the high-frequency shear mod- 
ulus A(lu A) does not vanish as it was assumed in 
that equation. Indeed, by adding a frequency indepen- 
dent term to Eq. (20) of Ref. we can recover all low-q 
predictions of Ref. [H without contradicting the stability 
criterion. 
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APPENDIX A: ELECTROSTATICS ON A 
CYLINDER 

In this appendix we calculate a few important ingredi- 
ents of a classical neutral jcllium model on a cylinder 

The interaction potential V(r — r') of two point unit 
charges living on a cylindrical surface is a periodic solu- 
tion of a 2D Poisson equation 



V 2 V(x, y) = -4nS(x - x' 
V(x + L,y) = V(x,y), 



nL)6(y-y') 



(Al) 



where L is the circumference of the cylinder and n is an 
integer number. This equation is readily solved by the 
Fourier transformation. The result takes the following 
form 

2 



V(r-r') = -2tt 



12/ — 2/ 



-In 



1 



- x' + i\y~y 



(A2) 
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It is easy to see that this potential has correct ID and 
2D asymptotic forms 



V(r 



-%\y-y'\, \y-y'\^L 

In Ir - r'l 2 , |r - r'l < L 



It is also straightforward to check the following useful 
representation for the potential V(r — r') 



t r / a y + y' , 

V(r — r ) = — 2tt — In 



^2iriz I 'L litiz /L 



(A3) 



Now we find the interaction potential U P b(r) of a nega- 
tively charged particle with a positive background charge 
density n b (r) that is homogeneously distributed on a 
cylinder and occupies the region —^P- < y < -^p- 



»Mr)=/>«(^-|y| 



where p — jj^r, and N is the total charge of the back- 
ground. The potential U pb (r) is given by the following 
integral 



U pb (r) = - J dr'V(r - r')n b {r') = -tt 



N Lt 2 
~2~T 



LT2 \ 



Lt 2 



- 2npy z + 2,p[\y\-^j 9 {\y\ - -f j (A4) 

Similarly one calculates the background-background 
interaction energy 



W bb = \J dvdr'V(r-v')n b (v')n b (v) = (A5) 
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